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The self-modulation instability is a key effect that makes possible the usage of nowadays proton 
beams as drivers for plasma wakefield acceleration. Development of the instability in uniform plas¬ 
mas and in plasmas with a small density up-step is numerically studied with the focus at nonlinear 
stages of beam evolution. The step parameters providing the strongest established wakefield are 
found, and the mechanism of stable bunch train formation is identified. 
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I. INTRODUCTION 

Plasma wakefield acceleration driven by proton beams 
is a concept that may open a path to future high energy 
lepton colliders [TJ. The interest to this concept is moti¬ 
vated by the ability of plasmas to support strong electric 
fields and the ability of proton beams to drive these fields 
over long distances. 

Available high energy proton beams consist of rather 
long 10 cm) bunches and cannot directly excite plasma 
waves of the required sub-millimeter length, as the beam 
frequency spectrum has a negligible component at the 
plasma frequency. Consequently, the proton bunches 
must be somehow adapted to the wave. Two meth¬ 
ods for that have now been proposed and studied. One 
method consists in longitudinal compression of the bunch 
as a whole, which is superior in terms of wakefield ex¬ 
citation EHs], but is hardly doable within the exist¬ 
ing infrastructure gHg. Another method involves beam 
micro-bunching in the plasma due to development of the 
self-modulation instability (SMI) [7], 10, |TT|. The lat¬ 
ter method is the cornerstone for the first experiment 
on proton driven plasma wakefield acceleration at CERN 
named AWAKE [U 1T2MT4] and for several supporting ex¬ 
periments with electron beams in other labs USHEj- 

The SMI appears due to mutual amplification of the 
beam radius rippling and the plasma wave [20] . The lin¬ 
ear phase of the wave growth has been studied in several 
papers and is already well understood. For the case of 
narrow (compared to the plasma skin depth), low density 
(compared to the plasma density), constant emittance 
and constant current beams, analytical expressions for 
the growth rate and phase velocity of the wave are avail¬ 
able Q3JI2IH23], which agree well with computer simula¬ 
tions. The SMI in slightly non-uniform plasmas is ana¬ 
lyzed in Refs. [2IJ [25] with the conclusion that the weak 
non-uniformity does not suppress the self-modulation, 
though the critical non-uniformity level differs by order of 
magnitude in theory [241 and simulations [25-. The insta¬ 
bility can grow either from an artificially introduced seed 
perturbation, or from the shot noise that is calculated in 
Ref. [26]. Various options are possible for the seed per¬ 
turbation: a short electron bunch [2H], a powerful laser 
pulse [27] . a sharp cut in the bunch current profile [18], 


or a relativistic ionization front co-propagating with the 
drive bunch [28]. Currently, the co-propagating relativis¬ 
tic front is considered as the mainstream variant, as this 
solution is chosen for AWAKE experiment. 

The seed perturbation becomes a must for beam slicing 
with the SMI. Otherwise several unstable modes could 
simultaneously grow and destroy the beam completely 
HU [20]. The most dangerous competing mode is the 
hosing (a non-axisymmetric mode belonging to the same 
instability family) in which the plasma wave couples with 
transverse displacements of beam slices. The interplay 
between hosing and SMI is studied in Refs. [294f3Tl and 
is by now limited to the analytically tractable case of 
narrow and constant emittance beams. No detailed nu¬ 
merical studies of mode interaction has been done yet, as 
this problem is essentially three-dimensional and compu¬ 
tationally heavy due to long beams and long interaction 
distances of interest. 

The nonlinear stage of SMI has been studied with nu¬ 
merical simulations only. It was shown that the wakefield 
amplitude cannot exceed the limit imposed by nonlinear 
elongation of the wave period [32] , and the designed op¬ 
eration regime of the AWAKE experiment is close to this 
limit [33] . A finite plasma radius or weak radial non¬ 
uniformity of the plasma have a small effect on the wave 
growth, if the radial scale of density variation is much 
greater than the plasma skin depth [34] [35]. As self- 
modulating beams are many plasma wavelengths long, 
motion of plasma ions could come into play and hamper 
wave excitation, if plasma ions are light [Ml El]. 

Once the self-modulation is completed, that is, the 
beam is split into well separated micro-bunches, the 
beam evolution continues towards some strongly de¬ 
stroyed state [38] . As a result, the wakefield amplitude as 
a function of the propagation distance drops down soon 
after reaching the maximum [39j. As a cure for beam de¬ 
struction, the small step in plasma density was proposed 
in Ref. [ 55 ] • With the density step, the wakefield ampli¬ 
tude remains at some high level for a long distance [39]. A 
proper longitudinal non-uniformity of the plasma density 
is thus a necessary element of a future plasma wakefield 
accelerator if the latter is based on self-modulation of 
long proton bunches. 

Optimum locations and magnitudes of density steps in 
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Refs. [33 39] were found by means of a straightforward 
numerical optimization without deep studies of how the 
step works and why its optimum is there. In this paper, 
we clarify mechanisms of beam destruction in the uniform 
plasma and the role played by the density step. Wherever 
possible, we compare two cases (the uniform plasma and 
the optimum density step) to emphasize the effect of the 
step on beam evolution. In SecjTIJwe define the problem 
under study. In Sec |III| we describe macroscopic proper¬ 
ties of beam evolution observed in simulations. Then we 
focus on the motion of beam particles in Sec |IV| and ana¬ 
lyze contributions of individual micro-bunches in SecfV| 
Main findings are summarized in Sec |VI| 


II. STATEMENT OF THE PROBLEM 


We use units of measure that are traditional for plasma 
wakefield studies: speed of light c for velocities, electron 
mass m for masses, initial plasma density no for den¬ 
sities, inverse plasma frequency uo~ x for times, plasma 
skin depth c/uj p for distances, and wavebreaking field 
Eq = mcuj p /e for fields. The SMI perturbs the beam 
axisymmetrically, so cylindrical coordinates (r, ip, z) and 
the co-moving coordinate £ = z — ct are convenient. 

As we aim at the detailed study of the nonlinear insta¬ 
bility stage, we intentionally exclude all irrelevant effects 
by the choice of beam and plasma parameters. In particu¬ 
lar, we consider the radially uniform cold plasma with im¬ 
mobile ions and the mono-energetic positron beam with 
the density profile 


n b0 e- r2/(2<T y £ < 0, 

0, £ > 0, 


( 1 ) 


where n b o = 4 x 10 -3 and ay = 0.5. The hard leading 
edge seeds the instability in the same way as in most 
discussed experiments. Unless stated otherwise, we sim¬ 
ulate about 25 wave periods, so the plasma response to 
the beam is surely linear for the chosen beam density. 
The positive charge of the beam makes discussions of 
potential wells more intuitive. The beam radius is cho¬ 
sen as a compromise between the analytically tractable 
case a r <C 1 and the design value of AWAKE experi¬ 
ment a r = 1. The beam relativistic factor 7 ^ = 1000 
is high enough to separate the betatron frequency and 
the plasma frequency. Using light particles (positrons) 
greatly speeds up simulations as compared to the proton 
beams. The beam angular spread Sa = 2 x 10 -4 is low 
enough to avoid the emittance-driven beam divergence 
that would otherwise complicate the self-modulation |39j . 

To ease comparison with earlier theoretical results, we 
also introduce the inverse betatron frequency 


Tb — y/Z'Yb/nbo ~ 707.1, (2) 

which is the unit of time, e.g., in Refs. [lTJ 21 ]. 

Since the plasma response is linear and axisymmetric, 
we simulate the SMI with the axisymmetric quasi-static 
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FIG. 1. On-axis electric field E Zjax and wakefield potential 
Tax of the plasma wave at two propagation distances: the 
beginning of the interaction at z = 0 (a) and a place of the 
almost maximum wakefield in the uniform plasma at z = 
1000 ~ y/2Tb (b). Regions of first five bunches are indicated 
by vertical lines and numbers. Definitions for and <£7 are 
illustrated by the example of the fifth bunch. 


2d3v hybrid code LCODE [40, [4l], in which the beam 
is modeled with macro-particles, and plasma is treated 
as the electron fluid. In some cases the result is cross¬ 
checked with the much slower fully kinetic version of 
LCODE [42j[43]. 

Typical patterns of the plasma wake are shown in 
Fig -0 where the subscripts “ax” denote on-axis values 
of the quantities. Along with the longitudinal electric 
field E z on the axis, we also show the on-axis value of 
the wakefield potential 

f‘OC> 

E z (r, ?)<%'. (3) 

The force acting on an axially moving ultrarelativistic 
particle is just (—VT). Regions of negative T are thus 
the potential wells that can confine beam particles. 

It is natural to consider parts of the beam which lie 
between adjacent maxima of T ax (0 as “bunches”. This 
definition is convenient in that the bunches can be iden¬ 
tified even at early stages of their formation. We define 
the bunch coordinate as the location of local poten¬ 
tial maximum at the “end” of the bunch and denote the 
wakefield amplitude after passage of the bunch by <£ 7 , 
where j is the bunch number. These definitions are il¬ 
lustrated in Fig.[ljb). The chosen definition of <£7 (more 
complex one than simply a value) is necessary for exclud¬ 
ing the slowly varying component of the wakefield which 
appears due to the average beam current. We also intro¬ 
duce the maximum wakefield amplitude $> m (z,L) that 
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FIG. 3. Growth of the wakefield amplitude with the beam 
length L and the number of passed bunches N at the moment 
of maximum wakefield (4> a ) and in the established equilibrium 
state ($/). 

that is typical for the SMI. 

In Fig. [2] we can distinguish five stages of beam evolu¬ 
tion in the uniform plasma. At the 1st stage, the wake¬ 
field structure changes from that of the seed perturbation 
to that corresponding to fastest instability growth. Then 
there is the 2nd stage characterized by the nearly expo¬ 
nential growth of the wakefield amplitude. This growth 
closely fits the theoretically predicted rate [21] 


FIG. 2. The maximum wakefield amplitude as a func¬ 
tion of propagation distance z for the uniform plasma (red 
lines) and the plasma with the optimum density step (blue 
lines) in normal (a) and enlarged semi-logarithmic (b) scales. 
The green dash line in (b) shows the theoretically predicted 
growth. The numbers and vertical lines indicate stages of 
beam evolution in the uniform plasma. Circles in (a) show 
points of absolute maximum <f> a and established amplitude 
$/■ 

we consider as a function of propagation distance z and 
beam length L. In what follows we put L = 160 unless a 
dependence on L is studied. 


III. PHENOMENOLOGY OF 
SELF-MODULATION 

In this section, we describe general SMI features and 
introduce the terminology necessary for further studies. 

First we consider the maximum wakefield amplitude 
<F m in relation to the propagation distance z (or the 
propagation time, which is the same) for the uniform 
plasma and the plasma with the optimum density step 
(Fig. [ 2 ]). Here the optimum step means a steep increase 
of the plasma density by 8.5% at z = 360. This com¬ 
bination of numbers maximizes the established wakefield 
amplitude that we measure at z = 10000. Quantitative 
characteristics of the presented graphs are beam param¬ 
eter dependent, so we focus on the qualitative behavior 


oc exp L 1/3 (A) j , (4) 

which is shown in Fig.[2|b) by the dash line. The 
2nd stage gives place to the 3rd stage (non-exponential 
growth) well before the maximum amplitude is reached. 
The 4th stage is that of the fast field decrease, and the 5th 
stage is characterized by the almost constant wakefield. 
With the density step, these stages are also observed, 
though the non-exponential growth is longer, and the es¬ 
tablished wakefield amplitude is higher. 

The graphs $> m (z) in Fig.[2ja) contain two significant 
points each: the absolute maximum of the amplitude 
<F a and the established amplitude <£>/. In the uniform 
plasma, both amplitudes monotonically grow with the 
beam length L (Fig. [ 3 ]) and, therefore, show how indi¬ 
vidual bunches participate in driving the wave. At the 
point of maximum field (<F a ), all available bunches work 
equally. At the established state (4>/), contributions of 
bunches decrease with the bunch number N. In other 
words, the longer the beam, the smaller the ratio of the 
established field to the peak field. There is no sense in 
plotting similar graphs for stepped density profiles, as 
the choice of the optimum case depends on the number 
of bunches. 

The purpose of the density step is to make the estab¬ 
lished wakefield amplitude as high as possible. We can 
vary location z s , magnitude Sn, and steepness of the step 
in search of maximum . For the sharp density step and 
L = 160, the function &f(Sn,z) is shown in Fig. [ 4 ] We 
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FIG. 6. The maximum wakefield amplitude <£> m as a function 
of propagation distance z (a) for selected longitudinal plasma 
profiles (b). 


FIG. 4. The map of step efficiencies $f(5n,z) in two projec- engineering formula for this dependence 
tions. 



FIG. 5. The magnitude of the optimum density step Sn as 
a function of the beam length L. Points are obtained from 
simulations, lines are approximations. 

see that the area of optimum parameters extends over a 
wide interval in covering the whole 2nd (exponential) 
stage of beam evolution [compare with Fig.[2][b)]. The 
wakefield amplitude reached before the step varies over 
an order of magnitude in this interval. The required mag¬ 
nitude of the step, however, is well defined for each step 
location. There are also several higher order maxima at 
multiples of the optimum step magnitude. 

The optimum magnitude of the plasma density step Sn 
depends on the beam length L or, equivalently, the num¬ 
ber of bunches N ~ L/(2tt) (Fig. [5]). There is a simple 


shown by the thick red line in Fig. [5] The inverse depen¬ 
dence on L is expectable, as it provides the same average 
shift of the wave with respect to beam for beams of dif¬ 
ferent lengths. The numerical factor in this dependence 
means that the optimum step makes the beam exactly 
one plasma period longer, if measured in local plasma 
wavelengths. The correction subtrahend in the denomi¬ 
nator indicates that the first two bunches play a special 
role in the self-modulation. Without this subtrahend, 
the approximation Sn = 2/N deviates from simulations 
at small L (thin green curve in Fig.[5|. Earlier simu¬ 
lations of low-emittance beams also fit the formula (J5|: 
the rightmost (cruciform) point is that for LHC beam 
j39| with cr^o taken as L. 

Figure [6](a) shows the time evolution of the maximum 
wakefield amplitude 4> m for several typical steps. These 
steps are drawn in scale in Fig.[6][b) and indicated in 
Fig.[4](b) with the same characters. We see that the max¬ 
imum field behaves similarly at different points of the 
“optimum plateau” and at higher order maxima. The 
smooth density increase (line “S”) produces almost the 
same result as the sharp step located at the beginning 
of the increase (line “B”). We therefore concentrate on 
steep steps, as a reasonably short length of the density 
transition has a little effect on the resulting field. 
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FIG. 7. Beam particles (points) and the separatrix (wireframe surface) in (r, £,p r ) space at different times for the uniform 
plasma case. In right fragments, small arrows show locations of the maximum on-axis beam density. 
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IV. MOTION OF BEAM PARTICLES 

The bunches are created and then partly destroyed be¬ 
cause of the radial motion of particles in potential wells 
of the wakefield. A slow longitudinal motion of poten¬ 
tial wells with respect to particles is also important. We 
therefore look at the motion of beam particles in the 
space (r, £,p r ), where p r is the radial momentum. The 
azimuthal angle cp is not important due to the axial sym¬ 
metry, the conserved azimuthal momentum p^ is negli¬ 
gible at low emittances, and the longitudinal momentum 
p z is almost constant at the considered time intervals. 
The potential wells in this space can be conveniently vi¬ 
sualized by the separatrix, which is the surface separating 
particles of positive and negative energy of the transverse 
motion. The separatrix is defined as 

p r = ±V'- 27 6 $(r,C at $(r,£)< 0 (6) 

and usually has the shape of a cone with negative-energy 
particles inside. We consider the interval —160 < £ < 
— 147.8 and display the resulting three-dimensional por¬ 
traits of the beam and the separatrix in two projections 
at several time slices. 

The case of the uniform plasma is shown in Fig.[7| At 
1st and 2nd stages of beam evolution [fragments (a)-(d)], 
the transverse motion of particles is field dominated. The 
separatrix width scales as the square root of the wave 
amplitude, while the particle transverse momentum is 
directly proportional to the amplitude, so the separatrix 
is much wider than the beam footprint at low amplitudes 
[Fig.[7^b,d)] . This means the transverse momentum of a 
particle is much smaller than that needed for escaping 
the potential well or than the momentum a particle can 
gain falling into the potential well. Since the wave am¬ 
plitude quickly grows during these stages, we can also 
neglect the momentum imparted to the particle at ear¬ 
lier times when the potential well had a different shape 
or location. Simply stated, at any time we can analyze 
the particle motion as if the particle recently appeared 
in the potential well with zero momentum. 

Regions of the decelerating field in the plasma wave 
are shifted a quarter-period with respect to potential 
wells. In the right column of Fig.[7| this is the right (for¬ 
ward) shift [see Fig.[7][b)]. The beam drives the wave only 
if there are more decelerated particles than accelerated 
ones; otherwise the wave would damp as the bunches pass 
by. Consequently, the potential wells must be populated 
asymmetrically with beam particles. The asymmetry is 
achieved by the backward motion of potential wells with 
respect to the beam (Fig.7j right column). The mecha¬ 
nism of the motion we will discuss later. With this mo¬ 
tion, the focusing force first tightens the beam in some 
place, and then the well falls back leaving the densest part 
of the beam (marked by small arrows in Fig. [7]) in the de¬ 
celerating field. This type of wave-particles interaction 
turns out to be the most efficient way of wave excitation, 
so it dominates over other possible perturbations and 
determines the overall beam evolution. The backward 
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FIG. 8. Bunch coordinates C as functions of the propagation 
distance z for the uniform plasma (red) and the optimum step 
(blue). Every second bunch is shown. 

motion of the wave was characterized in Refs. EUH2] in 
terms of the phase velocity. 

As the potential well moves with respect to the beam, 
it ingests beam particles at one side and releases at the 
other. This process continues while the transverse par¬ 
ticle motion is negligibly slow, but eventually comes to 
the end due to the increased particle momentum or ra¬ 
dial displacement of particles [Fig.[7|[e)-(h)]. Once this 
happens, the exponential growth (2nd stage) changes to 
the non-exponential growth (3rd stage). The reduction of 
particle inflow to the potential wells favors the asymme¬ 
try of well population, so the wave amplitude continues 
to grow. The growth rate, however, is smaller, as the 
unbalanced escape of particles from the well results in 
progressive well depopulation. 

The wakefield amplitude reaches its maximum when 
the decelerating phase of the wave approaches the parti¬ 
cle void area [Fig.|7|[i,j)]. At this moment, a considerable 
number of particles is in the strong decelerating field, but 
some of them are outside the separatrix. These particles 
gradually move aside from the axis causing the wakefield 
reduction during the 4th stage. This is a positive feed¬ 
back effect, as the reduced wakefield results in additional 
losses of particles from the well. Another reason for the 
extra particle loss is the motion of the potential well that 
continues during the 4th stage [Fig.[7|[j,l)]. At the stage 
of established wakefield (the 5th stage), potential wells 
and particle distributions inside the wells change a little 
compared to Fig. 0k, i). 

The motion of potential wells with respect to the co¬ 
moving frame is quantitatively characterized by Fig. [8] 
The motion with respect to beam particles is almost the 
same, as the average particle z- velocity (determined by 
the transverse momentum) differs from the light velocity 
by the term of the order of p*/( 2y^) < 10 -5 . In the 
uniform plasma, the motion stops at approximately the 
middle of the 4th (field decrease) stage. The integral shift 
is expectedly larger for tail bunches. 































7 



FIG. 9. Beam particles (points) and the separatrix (wireframe surface) in (r, £,p r ) space at different times for the optimum 
density step. In right fragments, small arrows show locations of the maximum on-axis beam density. 































FIG. 10. The effective beam current / e ff(0 at several propagation distances z in the uniform (a) and stepped-up (b) plasmas. 
The zero current is at the bottom of each drawing zone. Blue parts of the curves correspond to focusing phases of the wakefield 
($ ax < 0), orange parts to defocusing (<h ax > 0). Horizontal dotted lines show the effective current of the initial beam. Red 
vertical lines help to relate locations of bunches at different times with their final locations. 


Beam and separatrix portraits for the case of the op¬ 
timum density step (profile “A” in Fig. [6]) are shown in 
Fig. § Snapshots (a) and (b) are taken right before the 
density step and are the same as in Fig.[7|(c,d) , fragments 
(c) and (d) are taken soon after. The step makes the 
plasma wavelength shorter, but does not change the num¬ 
ber of wakefield periods (Fig. [8]). The partially bunched 
(or rather, rippled) beam thus drives the wave with the 
period that is substantially longer than the local plasma 
wave period. The increased density acts as a “spring” re¬ 
sisting further elongation of the wakefield period. If the 
“spring pressure” is properly adjusted by the step magni¬ 
tude, then the wake period does not change, and poten¬ 
tial wells do not move with respect to the beam (Fig.|8j 
right column of Fig.[9|. The strength of the “spring”, 
however, depends on the wave amplitude: the higher the 
stronger. This is because the own wakefield of a particle 
bunch has a little effect on the high amplitude wave, but 
easily modifies a weak wave. We therefore observe some 
backward shift of the wave in the time interval of the 
low-amplitude wave (z < 800 in Fig.|8|. 

There is, however, the question of how to provide the 
asymmetry of well population necessary for efficient driv¬ 
ing the wave. As we see from Fig.[9j the backward shift of 
the wave at z ~ 600 is responsible for that. The shifted 


potential well imparts a large transverse momentum to 
some beam particles [Fig.|9|e,f)] that fall outside the sep¬ 
aratrix when the well shifts forward [Fig.[9^g-j)] and cre¬ 
ate a gap in the established well population [Fig.[9jk,l)]. 

Evolution and phasing of the bunch train as a whole 
can be examined with the help of the effective current 

rOO 

l e s(0 = / n b (r , £)e _r 2irr dr (7) 

Jo 

that qualitatively accounts for reduced contribution of 
defocused particles to the wakefield. Strictly speaking, 
the contribution of an off-axis particle to the on-axis 
electric field E z is proportional to the modified Bessel 
function of the second kind K${r) [26} EH- However, the 
exact formula is impractical when applied to simulations 
output because of the noise generated by near-axis parti¬ 
cles, so we use 0. it is of particular interest to correlate 
the effective current and the location of focusing and de¬ 
focusing regions. 

The map of effective beam currents for the uniform 
plasma [Fig.[lQ|^a)] reveals that all bunches in the beam, 
except the first one, behave similarly: leading edges of 
the bunches are in defocusing phases, trailing edges are 
focused, and the bunches move backward in £ until their 
trailing edges run into zero I e q regions. After that, the 
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defocusing phases “eat up” large parts of the bunches, 
the larger the further the bunch is from the beam head. 

In the stepped-up plasma [Fig.[TQ^b)] ? we observe the 
transition between two typical cases as we go from head 
to tail of the beam. For the head bunches, the con¬ 
sidered density step is smaller than the optimum one, 
so the bunches evolve almost as in the uniform plasma: 
the final bunch shapes at £ > —40 in fragments (a) and 
(b) of Fig. [lO] are similar, and the bunches move back¬ 
ward at the nearly the same speed. For the tail bunches 
(£ ^ —160), the backward motion is suppressed by the 
reduced plasma wavelength, but the bunches in the final 
state fill potential wells almost symmetrically and drive 
the wave inefficiently. In the middle (£ ^ —100), the final 
bunch shape looks optimal, as decelerated particles sur¬ 
vive, while accelerating phases of the wave are depleted. 

As we see, the density step instantly switches the 
beam-plasma system from the fastest growing mode to 
the state characterized by a slower wave growth and a 
favorable final beam state. This is possible owing to the 
ability of the grown SMI mode to suppress development 
of other modes. The same feature enables suppression of 
the hosing in the case of a seeded SMI. Otherwise, growth 
of the detuned mode would be hampered by the fastest 
mode at the new plasma density. 

The optimum switching occurs at the linear stage of 
the wave growth at which the ratio of wave amplitudes 
before and after the step does not depend on the am¬ 
plitudes themselves. This is the reason why the interval 
of possible step locations is as long as the whole second 
stage of beam evolution. 


V. CONTRIBUTIONS OF SEPARATE 
BUNCHES 


In the previous section, we looked at the process from 
the beam’s standpoint and mostly examined how the 
beam is changed by the wakefield. Here we take the 
plasma wave as a reference and consider the complemen¬ 
tary part of the process, the wakefield excitation by the 
bunched beam. 

We start with reminding few important consequences 
of the linear wakefield theory [45]. A short particle bunch 
traveling in the plasma wave modifies this wave by its 
own wakefield (Fig. 11). The result of modification de¬ 
pends on which phase of the wave the bunch is located 
at. If the bunch is at the maximum decelerating gradi¬ 
ent, its wakefield is in phase with the master wave, and 
the bunch simply amplifies the wave [Fig.JlTJb)]. If the 
bunch is located at the bottom of the potential well, the 
wave amplitude remains unchanged, but the phase is ad¬ 
vanced [Fig.[TT]^c)]. At intermediate locations, the bunch 
simultaneously amplifies the wave and moves it backward 
[Fig®)]- Note that same conclusions are true for elon¬ 
gated symmetric bunches like that shown by the dash line 
in Fig.[Hp). 

The above considerations explain why the wave motion 



FIG. 11. Illustration of wakefield modification by a short 
bunch. Wakefield of the bunch is: (a) drawn separately, (b) in 
phase with the master wave, (c) approximately 7 t/ 2 delayed, 
and (d) 7 t/4 delayed, (e) Composition of complex amplitudes 
for the three cases. 


with respect to the beam continues after reaching the 
maximum wave amplitude in the uniform plasma. At the 
maximum amplitude, front parts of the bunches are in the 
defocusing field [Fig. T^i,j), Fig. 10 [a)]. These parts must 
disappear and thereby shift bunch centroids backward, 
thus shifting the wave backward. 

The wave modification by individual bunches can be 
conveniently characterized by the complex amplitude 


Fj = $ i exp(-i^V«), 


(8) 


where j is the bunch number, and n is the local plasma 
density that affects the local plasma wavelength. Each 
bunch makes the increment A Fj = Fj — Fj- 1 to the 
complex amplitude. It is more informative, however, to 
look at the quantity 


dFj = A Fj exp(i£j_i y/n) (9) 


that characterizes the bunch contribution in relation to 
the wave excited by previous bunches [Fig.JTlJe)]. 

If the bunch is a fragment of some constant-current 
long beam, then its contribution to the existing wakefield 
is strongest if the bunch completely fills the decelerating 
and focusing phase of the wave (Fig.[l2|). A shorter or 
partially destroyed bunch is obviously less efficient. A 



fewer useful particles optimum 


FIG. 12. Illustration of the optimum bunch positioning. 
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FIG. 13. The complex amplitude Fj (a,d) and the increments dFj to the complex amplitude (b,c,e,f) at different propagation 
distances for uniform (a,b,c) and stepped-up (d,e,f) plasmas. Points of the same color correspond to different bunch numbers j, 
lines between them are drawn for better visibility. In fragments (b,c,e,f), small circles mark contributions of the second bunch, 
points on the abscissa axis are contributions of the first bunch, gray shadings show areas of defocusing risk, and dashed lines 
are theoretical approximations explained in the text. 


longer bunch contains either accelerated particles that 
reduce the wave or defocused particles that make the 
bunch unstable. The optimum angle between dFj and 
Fj _i is thus 7 t/4 , as in Fig.: 
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d). 


Growth of complex amplitudes at different propaga¬ 
tion distances is shown in Fig.[l3ja,d) for the two studied 
cases. With the density step, the spirals are expectedly 
tighter, since the ratio of the bunch-to-bunch distance to 
the plasma period is higher in this case. Other differences 
are seen only at the graphs of individual bunch contri¬ 
butions. In the uniform plasma, the five stages of beam 
evolution correspond to qualitatively different behavior 
of dFj. At the first stage (z < 200), the bunches appear 
where the potential wells of the seed wakefield are, and 
their contributions are purely imaginary: no amplitude 
increase, only the phase shift [Fig.[l3|[b)]. The exponen¬ 
tial growth stage (200 < z < 600) is characterized by 
the nearly constant phase shift (imaginary part) and the 
amplitude increment (real part of dFj ) that grows along 
the beam. The transition from the exponential to the 
non-exponential growths (z ~ 600) occurs when contri¬ 
butions of tail bunches approach that of a fully bunched 
beam (\dFj\ ~ 0.002 in our case). As soon as the bunch 
contribution cannot grow further in absolute value, it 
moves down to the area of defocusing danger [gray areas 
in Fig. 13[b,c)]. In other words, the bunch advances the 
wave phase less and therefore increases the local phase 
velocity of the wave. Staying in the gray area is accom¬ 


panied with loss of particles, so the bunch contribution 
reduces in the absolute value at z > 1000 [Fig.|l3|[c)] 
until dFj exits the defocusing area. The established con¬ 
tributions of the bunches (at 2 = 10000) are thus optimal 
in phase, but small in absolute value, so the established 
wakefield is relatively weak. 

The picture of bunch contributions after the density 
step is less clear and even looks chaotic [Fig.[l3](e,f )]. The 
general trend is that all contributions are oriented to¬ 
wards the vertical axis, and the main job of the bunches is 
to push the wave backward. Some of the bunches are ac¬ 
celerated [reside in left half-planes of Fig.[l3|[e,f)], but the 
decelerated ones are more numerous, so the wave grows, 
in average. Importantly, bunch contributions never en¬ 
ter the defocusing area and remain strong in absolute 
value. It is also worth mention that, in the established 
state [z = 10000 in Fig.p^f)], the last bunches are not 
the most efficient ones, as we noticed when discussing 
Fig. [lO] The general conclusion from this picture is that 
bunch contributions in the case of the optimized density 
step are far from being regular and optimum, and there is 
a room for improvements with more sophisticated density 
profiles. 

Let us relate the observed bunch contributions with 
theoretical models. There are two of them available. One 
model is the wave growth in the linearly responding uni¬ 
form plasma m- It is characterized by the amplitude 
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(a) 




FIG. 14. The complex amplitude Fj and its theoretical approximations for the uniform plasma at z — 500 (a) and z — 10000 
(b) and for the optimum density step at z — 10000 (c). 



FIG. 15. The established complex amplitude Fj at z = 10000 
for the uniform plasma and several density steps labeled ac¬ 
cording to Fig. [6] 


growth @ and the phase shift 


Scpj = -£j - 2irj = 


3 L 1/3 


nj 


2/3 


( 10 ) 


which can be combined to yield the spiral 


oc exp (a Stfij) , Fj oc exp ((a + i)5(/)j) (11) 

with a = y/3 being the slope of the spiral with respect to 
a constant radius circle. Another model [46] is the opti¬ 
mum train of rectangular bunches that completely fills fo¬ 
cusing and decelerating phases of the wave. In the limit of 
large bunch numbers it forms the spiral (11) with a = 1. 


The theoretical spirals are close to numerically obtained 
ones at corresponding times (Fig.[l4|), though more sen¬ 
sitive diagnostics (maps of dFj) reveals deviations from 
expected slopes [dashed lines in Fig.[l3|[b,c)]. The estab¬ 
lished wake in the case of the density step also follows 
the spiral but with a ~ 0.6 [Fig.[f4|[c), Fig.[l3|[e)]. 

At multiples of the optimum step magnitude, some 
bunches are almost completely destroyed and make no 
contribution to the wakefield, while others resonantly 
drive a substantial wake (Fig.[l5]). This effect is simi¬ 
lar to wakefield excitation by a non-resonant bunch train 
m and involves partial or complete defocusing of the 
bunches which contribution to the wake is destructive. 


VI. SUMMARY 


Let us summarize the key findings of the paper. 


In the uniform plasma, the seeded self-modulation in¬ 
stability first transforms the long particle beam into the 
train of well separated bunches and then partly destroys 
these bunches. The larger the bunch number, the smaller 
part of the bunch remains. Both creation and destruction 
of the bunches occur as a result of particle defocusing by 
the plasma wakefield. Development of the SMI is accom¬ 
panied with a flow of beam particles through potential 
wells of the plasma wave. The wells crawl backward along 
the beam, grab new portions of beam particles from one 
side, and release particles from another. As soon as the 
potential wells run into particle-depleted cross-sections 
of the beam, the SMI passes into the nonlinear stage, 
and the wakefield growth slows down. At the propaga¬ 
tion distance where the maximum wakefield is reached, 
the bunches are well separated in space and coherently 
drive the wave. The wave motion with respect to the 
beam, however, continues after this point and defocus 
major parts of the bunches. The remains of the bunches 
are optimally spaced for stable wakefield excitation, but 
are too small to excite a strong wave. 


The character of beam evolution can be changed by 
increasing the plasma density at the linear stage of insta¬ 
bility growth. The optimum magnitude of the increase is 
such that to make the beam one wavelength longer if mea¬ 
sured in local plasma wavelengths. In the plasma of the 
increased density, the bunch-to-bunch distance is longer 
than the plasma wavelength, and the average phase veloc¬ 
ity of the wave equals the beam velocity. Nevertheless, 
well separated bunches are formed, and these bunches 
quickly come to equilibrium with the wave and excite 
the strong wakefield over a long propagation distance. 
The detailed study of the equilibrium bunch train, how¬ 
ever, reveals that the wakefield is excited not optimally, 
so there is a hope that even stronger waves can be excited 
with more sophisticated plasma density profiles. 
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